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A number of rotation-powered millisecond pulsars are powerful sources of X-ray 
emission. We present predictions for the spectral characteristics of these sources 
at gamma-ray energies, using a model for acceleration and pair cascades on open 
magnetic field lines above the polar caps. Since these pulsars have low surface mag- 
netic fields, the majority do not produce sufficient pairs to completely screen the 
accelerating electric field allowing particle acceleration to high altitude. The result- 
ing emission above 1 GeV comes from curvature radiation by primary electrons with 
radiation-reaction-limited Lorentz factors. The spectra are very hard power-laws with 
exponential cutoffs between 1 and 50 GeV, and the spectral power peaks near the cut- 
off energy. Millisecond pulsars are thus ideal targets for air-Cherenkov detectors that 
are able to reach energy thresholds below 50 GeV. 



1 Introduction 

There are more than 100 radio pulsars now known with periods between 1 and 30 
milliseconds. They are believed to have been spun-up by accretion torques due to 
mass flow from a binary companion 1 1 1 and have surface magnetic fields, in the range 
of 10 8 — 10 10 G, much lower than those of the normal pulsar population. Most are 
found in binary systems with low mass or compact (White Dwarf or neutron star) 
companions. Although they are found throughout the Galaxy, globular clusters seem 
to be particularly good breeding grounds for millisecond pulsars (MSPs) and a num- 
ber of clusters contain as many as 15 (in 47 Tuc) or 21 (in Ter 5) 1131 . Many MSPs 
are also seen at high energy wavelengths. About 30 MSPs are detected as X-ray point 
sources and pulsed emission has been detected in 7 of these (2||8]. Pulsed 7-ray emis- 
sion has been detected from one MSP, PSR J02 18+4232, by EGRET on the Compton 
Gamma-Ray Observatory 1 9 1 . Thus, MSPs must be powerful particle accelerators and 
despite their low surface magnetic fields, their spin-down luminosities are comparable 
to those of young pulsars. 

Polar cap particle acceleration theory (3] Q75| [5] |6 1 predicts that MSPs should 
accelerate particles to high Lorentz factors (~ 10 7 ). These high energies, combined 
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with the small radius of curvature of the open field lines, enables curvature radiation 
emitted by these particles to reach 10 - 100 GeV, depending on the pulsar parameters. 
Most should be easily detectable by GLAST and some may be detectable by Air 
Cherenkov Telescopes (ACTs). We will discuss how the predicted spectrum and flux 
depend on the period, magnetic field, inclination and equation of state of the pulsar 
and the prospects for detection of the emission by various instruments. 



2 Curvature Radiation from Accelerated Particles 

MSPs have the lowest magnetic fields of any radio pulsars and such low fields have 
two important consequences for their high-energy emission. They produce relatively 
few electron-positron pairs by one-photon pair production, a process that requires 
fields nearer the critical value of B cr = 4.4 x 10 13 G. As a result, the pair cascades 
in most MSP magnetospheres will not screen the electric field as they do in younger 
pulsars. Also, the spectral cutoff due to pair production attenuation, which depends 
inversely on field strength, will occur at much higher energy. Therefore, particle 
energies will be limited by curvature radiation reaction to a value 1 10][5]|6) 

where Eu is the parallel electric field and p c is the field line radius of curvature. 
The parallel electric field, including the contribution from inertial frame-dragging, is 
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where B% = B/10 8 G is surface magnetic field strength, P ms = P/l ms is pulsar 
period, x is the inclination angle, <j> is magnetic azimuth angle, 9q = {2-kR/cP) 1 ' 2 is 
the polar cap half-angle, rj = r/R is the radial coordinate and «o.i5 = O.I5/45/-R5, 
J45 = J/10 45 g cm 2 and R§ = R/10 6 cm are the moment or inertia and radius of the 
neutron star. In this case, the curvature radiation peak energy, SpJ^ k = (4/3)ecfl will 
be 



where Ac is the electron Compton wavelength. Since the polar cap half-angle 9q is 
large for MSPs, the radius of curvature of open field lines, p c ~ r/0o, is relatively 
small. Thus the CR spectrum of MSPs can extended beyond 10 GeV, but since the 
pair creation will attenuate the emission above e CS c ~ 10 6 Pn&Bg 1 , where e csc is 
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Figure 1 : Predicted curvature radiation (CR) and unattenuated (light line) and attenuated (dark line) inverse 
Compton (ICS) spectra of PSR B 1 937+2 1 , for two different inclination angles \ m radians and FP equation 

of state (see Fig. 2). 



the photon escape energy in mc 2 , the maximum energy is effectively limited to 100 
GeV or below. 

It should be noted that any detected emission above 2 GeV from MSPs cannot be 
produced in outer gap accelerators 1161 171. since the quantity E^^p]/ 2 determining 
the peak energy is smaller. Therefore, emission detected from MSPs at very high 
energies must originate from near the polar cap. 



3 Predicted Spectrum 

The predicted curvature radiation spectrum of electrons with radiation-reaction lim- 
ited Lorentz factors is very hard, with photon power law index —2/3 and exponential 
cutoff cxp (—e/ecR)- Figure 1 shows the calculated spectrum of curvature radiation 
for PSR B1937+21, which is the fastest MSP with period of 1.6 ms. It has a pre- 
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Figure 2: Predicted curvature radiation (CR) and unattenuated (light line) and attenuated (dark line) in- 
verse Compton (ICS) spectra of PSR J0437-4715, for soft (BPS), medium (FP) and hard (PS) neutron star 
equations of state (see [ 5 ]) and inclination angle \ = 0.1. EGRET upper limits are from (B). 

dieted CR peak energy around 25 GeV, which is not attenuated by pair creation since 
£csc ~ 30 GeV, so that significant emission extends above 50 GeV. The spectrum is 
shown for two different inclination angles and since the magnetic axis of this pulsar 
may be highly inclined, from the fact that both the radio and X-ray profiles have in- 
terpeaks, it may be more difficult to detect with H.E.S.S. or VERITAS. A component 
due to inverse Compton scattering (ICS) of thermal X-rays from the neutron star sur- 
face by accelerated electrons is also expected and will peak at energies of 1-10 TeV 
However, the ICS emission will be strongly attenuated by magnetic pair production 
to a level that is far below present detector sensitivities. 

Figure 2 shows the calculated spectrum of CR for PSR J0437-4715, illustrating 
the dependence of the emission on the neutron star equation of state (EOS). The values 
of i?6 and J45 for the different EOS determine the value of the k parameter in the 
parallel electric field of eqn 0. The apparent presence of a core component in the 



4 



B1 937+21 
B1821-24 _ 
B1 957+20 _ 
J021 8+4232 _ 

B1 257+1 2 
J1843-1 113 
J 1909-3744 
J0034-0534 
J0737-3039A 
J21 24-3358 
J 1024-0719 _ 
J21 29-5721 
J191 1-1114 
J0437-471 5 



10 



MAGIC 



H.E.S.S. 
VERITAS 



10 u 



10 1 



10* 



10 d 



Flux(>50GeV) (10 - 12 cnrf 2 s" 1 ) 



10 4 



Figure 3: Predicted range of integral flux (> 50 GeV) from individual millisecond pulsars. Light gray 
vertical bands are integral flux thresholds at > 50 GeV for H.E.S.S., VERITAS and MAGIC II 



radio profile suggests a small viewing impact angle to the magnetic pole. Therefore, if 
the acceleration model described above is correct, we should see the full low-altitude 
CR component. The predicted spectrum for the stiffer EOS (PS) is not consistent with 
the EGRET upper limits, suggesting that either the EOS is soft or that the model is 
incorrect. In the former case, PSR J0437-4715 will not be detectable by H.E.S.S. 



4 Predicted Flux From MSPs in the Galactic Plane 

Figure 3 displays a rank-ordered list of MSPs and their predicted integral fluxes > 50 
GeV, where the flux range reflects the dependence on inclination angle. Only a few 
MSPs, including B 1937+21, B 182 1-24 and B 1957+20, have a flux as predicted by 
this model that would be detectable above 50 GeV and all three have pulsed X-ray 
emission. Actually, B 182 1-24 is the only known MSP in the globular cluster M28, 
but at a distance of 4.9 kpc may be individually detectable. A few more MSPs would 
be detectable by ACTs such as MAGIC that may reach energy thresholds as low 
as 30 GeV and that may have lower flux sensitivity above 50 GeV. J02 18+4232 is 
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Figure 4: Predicted integral flux (> 50 GeV) from globular clusters containing millisecond pulsars. The 
number of known radio pulsars are shown in parentheses for each cluster. 



the only MSP detected above 100 MeV by EGRET |9| and may be detectable by 
ACTs depending on the viewing angle. J0737-3039A is the MSP with P = 23 ms 
in the double pulsar binary system 1 1 1 1 and whose radio emission is eclipsed by the 
magnetosphere of its 2.8 s pulsar companion. 



5 Predicted Flux From Globular Clusters 

Although globular clusters are generally more distant than many individual Galactic 
MSPs, the combined flux from multiple MSPs in these clusters may be detectable 
by ACTs. Figure 4 shows several globular clusters containing one or more MSPs, 
rank-ordered by their predicted combined flux > 50 GeV from the known radio pul- 
sars in the cluster. In particular 47 Tuc, containing fifteen known radio MSPs and 
possibly many more undiscovered MSPs, may be detectable by ACTs. The clusters 
NGC6397 and NGC6624 each contain only one known MSP, J1740-5340 and J1823- 
3021 A respectively, that individually may produce a detectable flux. However, there 
is a greater degree of uncertainty in predicting 7-ray flux from MSPs in globular clus- 
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ters since the measured period derivatives have model-dependent contributions from 
acceleration in the cluster potential 1141 . 

6 Conclusions 

Pulsed emission from several millisecond pulsars may be detectable by ACTs at en- 
ergy greater than 50 GeV, if emission from a low-altitude accelerator is visible. In 
particular, telescopes that are able to reach thresholds below 50 GeV should be able 
to detect all the X-ray bright MSPs. For large viewing angle to the magnetic axis, the 
low-altitude emission which reaches the highest energies will not be visible. In this 
case, even if emission from an outer gap accelerator is visible, no emission above 
10 GeV is expected. Several globular clusters, including 47 Tuc, NGC6397 and 
NGC6266, may also be detectable by ACTs. Detection or limits for individual pulsars 
believed to have small viewing angles to the magnetic pole, such as B 1937+21 and 
J0437-4715 will give direct constraints on the electrodynamics of particle accelera- 
tion. 
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